Local density of states in a d-wave superconductor with stripe-like modulations and a 

strong impurity 
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Using an effective Hamiltonian with d-wave superconductivity (dSC) and competing antiferro- 
magnetic (AF) interactions, we show that weak and one-dimensionally modulated dSC, spin density 
wave (SDW) and charge density wave (CDW) could coexist in the ground state configuration. With 
proper parameters, the SDW order exhibits a period of 8a, while for dSC and CDW orders the period 
is 4a. The local density of states (LDOS), which probing the behavior of quasiparticle excitations, 
is found to have the identical stripe-like structure as those in dSC and CDW orders. The LDOS as 
a function of the bias voltage are showing two small bumps within the superconducting coherence 
peaks, a signature of the presence of stripes. When a strong impurity like Zn is placed in such a 
system, the LDOS at its nearest neighboring sites are suppressed at the zero bias by the local AF 
order and show a double-peak structure. 

PACS numbers: 



Many of the anomalous properties of high-T c supercon- 
ductors (HTS) are believed to be related to the compe- 
tition between the <i-wave superconductivity (dSC) and 
the hidden antiferromagnetic (AF) order, particularly 
in the underdoped HTS . Inelastic neutron scattering 
(INS) studies of the magnetic fluctuations in some of the 
HTS samples have provided important clues to the na- 
ture of the electronic correlations within the doped CuC>2 
planesiS^. In additon, the existence of one-dimensional 
charge density wave (CDW) and spin density wave 
(SDW) has also been reported on Lai.e-xNdo^Sr^CuO^ 
and YBa 2 Cu 3 6 .3^. 

The elastic neutron scattering measurements on 
La2_ a; Sr a ;Cu04 with x = 0.10 sample by Lake et al& 
found that the signal got enhanced at (i, | + <5)(— ) at 
low temperature, indicating the presence of static AF or- 
der at least in some parts of the sample. The coexistence 
of the CDW and dSC in HTS in the absence of a magnetic 
field has attracted a lot of theoretical works2&2iiSiii. Re- 
cent progress in scanning tunneling microscopy (STM) 
on the surface of Bi 2 Sr2CaCu208+5 (BSCCO) has given 
us a high resolution probe of its electronic correlations. 
The STM imaging by Hoffman et ali— in the presence of 
a magnetic field revealed the quasiparticle states around 
the vortex cores in slightly overdoped BSCCO to exhibit 
a CDW of "checkerboard" pattern with 4a periodicity. 
Howald et alM^ reported that the four unit-lattice pe- 
riodic charge modulation survives even at zero magnetic 
field and is energy independent. And it has been at- 
tributed to the superposition of stripe phase a 13 ' 15 ori- 
ented along x- and y- directions in the sample. This 
issue has also been examined by several theoretical arti- 
cles in the presence of a magnetic fieldiLiSiiLiSiiii and in 
the absence of a magnetic fieldii. 

In case a stripe phase indeed exists in some of the 
HTS samples close to optimal doping as reported in the 
experimentaiii4, it should have important impact on the 
physics of HTS. We need to understand the origin of the 
stripe-like modulation observed in the local density of 
states (LDOS) and to examine its consequence on other 



experimentally measurable quantities, such as the LDOS 
near a unitary impurity in the absence of a magnetic 
field. In view of neutron scattering experiments men- 
tioned above, we assume that the stripe phase observed 
in the STM experiments is originated in the competing 
AF interaction. In this paper and based upon a model 
Hamiltonian, we first construct a superconducting phase 
for a nearly optimal-doped or slightly underdocd HTS 
sample in which the dSC coexists with the SDW, and 
CDW orders. The parameters will be chosen in such 
a way that the dSC and CDW will have a weak onc- 
dimension-like modulation with period 4a while the pe- 
riod for the SDW is 8a. The superposition of the x- and 
y- oriented CDW stripes should yield the experimentally 
observed checkerboard patterns! 3 ^. Then we calculate 
the LDOS of the stripes and compare it with the pure d- 
wave case. On the other hand, the low temperature STM 
experiment™ observed a sharp resonance peak near the 
zero bias at a nonmagnetic unitary impurity site, consis- 
tent with prediction^! for a pure dSC. In the presence 
of the stripe phase, we found that the LDOS is strongly 
suppressed by the AF order at the zero bias and a double- 
peak shows up. Hopefully, the obtained result could be 
used to compare with future STM experiments for sam- 
ples with magnetic originated stripes. 

Based on a square lattice with lattice constant 
a(=l),we start from the effective mean-field Hamiltonian 
with g as the nearest neighbor attractive d-wave pair- 
ing potential, and U as the on-site Coulomb interaction 
representing the competing AF order. 



H = -E*« c i<T c j- + E( A « c iT c ji + A y c ji c iT) 

ij.o- ij 

+X^( m ^-M+vocL c i* (i) 



where ty is the hopping integral, Ay = f ( c iT c ji — c U c jt) 
is the spin-singlet d-wave bond order parameter, ra- la = 
U{riia) is the AF order parameter, jj, is the chemical 
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potential. In addition there is also a nonmagnetic im- 
purity with V\ = VqSqi as the single-site scattering po- 
tential at site (0,0). We shall diagonalize the above 
Hamiltonian by using Bogoliubov transformation, c- la = 

X^[ M kr7no- — crv icrln&\ > an d the equations of motion 
for c\ a and q- will lead to usual Bogoliubov-de Gennes' 
equations (BdG), 
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where "Hij<y = — ^i+e,j + (mi ff — p)5i- } + Vb^oj- The sub- 
script e are the vectors x, y and x + y toward the nearest 
neighbor (NN) and next-nearest neighbor (NNN) sites, 
respectively. To self-consistently solve BdG equations we 
could get the N positive eigenvalues E n (n = 1 • • • N) 
and the N negative eigenvalues E n with corresponding 
eigenvectors {u^,v^) and (— , u™j*), respectively. Us- 



ing the following convenient notation = (- 



and v-\ = (u™?,vV\), the self-consistent conditions be- 



come, 
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A ij = El(«»+«)^nh(^) 



(3) 



where f(E) = l/(e /3E + 1) is Fermi distribution function. 

Since the calculation will be performed near the 
optimal-doped or slightly undcrdoped regime, we choose 

the filling factor, which is defined as rif = X^i<r( c L Cicr )/ 
N x N y with the summation over one unit cell, is fixed to 
be 0.85 (i.e. the hole doping Xh = 0.15), where N x , N y 
are the linear dimension of the unit cell under consider- 
ation. The chemical potential p, therefore, needs to be 
adjusted each time when the on-site repulsion U is var- 
ied. The NNN hopping integral is chosen to be t' = —0.2, 
as relevant to the hole-doped cuprate, to fit the hole-like 
Fermi surface. Through out this paper, we are going to 
use the same value of U, g, t', and rif. 

Once the self-consistent solution is obtained, we calcu- 
late the staggered magnetization and the electron density 
defined as M s = (— l) 1 ^^ — n-^) and rii = (n-^ + n-^), 
respectively. Furthermore, the LDOS of the energy E at 
the positon i can be written as 
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where Pi{E) is proportional to the local differential tun- 
neling conductance as measured by STM experiment, and 
the summation is averaged over a M x x M y wavevectors 
in first Brillouin Zone. In addition, the LDOS spatial 




FIG. 1: The spatial distribution of (a) The dSC order pa- 
rameter 5|Ai| = |Ai| — 0.041, (b) The staggered magneti- 
zation ,and (c) the electron density &n\ = J^o- 71 " " — n f f° r 
the pure superconductor (without impurity), (d) The LDOS 
maps Spi{E) = pi(E) -0.084 at the energy E = 0.04. The size 
of the unit cell is N x x N y = 32 x 32. The attractive potential 
is g = 0.71. The on-site repulsion is U = 2.05. The NNN 
hopping term is t' = —0.2. The filling factor is rif = 0.85. 



maps observed in STM experiment at a fixed bias volt- 
age E could be obtained by calculating pi (E) at each site 
of the lattice. 

In the absence of the magnetic field, ty = t and t' are 
the NN and NNN hopping integrals, respectively. Since 
the strength of the onsite repulsion, U , required to gener- 
ate AF order has to be fine tuned to study the interplay 
of the two competing orders. The non-zero AF order ap- 
pears only for certain range of U and g, and larger U/g is 
able to produce stronger antiferromagnetism. With our 
chosen parameters U = 2.05 and g = 0.71, weak stripe 
modulations can be generated in the dSC order param- 
eter, the staggered magnetization and electron density. 
Here we set t = 1. 

In Figure 1(a), our numerical result for the spatial dis- 
tribution of the dSC order, which is defined as Ai = 
[Aj + £ + Aj_£ — Ai + .y — Ai_g]/4 , exhibits the stripe be- 
havior along y-direction with the period 4a. The stag- 
gered magnetization (Fig. 1(b)) shows stripe like SDW 
along y-direction. Its period is 8a and its amplitude is 
less than 0.1. The electron density (Fig. 1(c)) has one- 
dimensional CDW modulation with 4a as its period and 
a very weak amplitude (less than 0.001). The origin of 
such static stripes could be understood in terms of the 
existence of a nesting wave vector qA ~ 0.257r/a connect- 
ing the upper and lower pieces of the Fermi surface near 
( 4 ) (tt,0) along fcy-direction [See Fig. 1 in Refs. 24]. For 
proper values of U and doping, this wave vector would 
modulate the staggered magnetization with SDW stripes 
along x-direction with period I'njqA = 8a. Accompa- 
nying the SDW stripes we have the charge stripes with 
period 4a. If we choose a larger U(> 2.05) here, the effect 
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is to enhance the amplitudes of the SDW/CDW stripes, 
not to change its period. For smaller U(= 2.0), only dSC 
without stripes is obtained. The stripe configurations 
discussed above are associated with the ground state of 
the system. What have been observed by STM exper- 
iments are related to quasiparticle excitations. In Fig. 
1(d) we present the spatial map of the LDOS at energy 
E = 0.04, and it exhibits the same stripe like modulation 
with period 4a as that appeared in CDW shown in Fig 
1 (c) . We have checked the LDOS at several different bias 
energies, and found that the same stripe like structure 
still prevails and the modulation period or wavevector is 
energy independent. This result can be regarded as a 
necessary condition for existing a static-stripe CDW or- 
der in the sample. Since the y- and x- oriented stripes are 
degenerated in energy. The observed checkerboard pat- 
tern in the STM experimentsi&ii could be understood 
by a superposition of these two perpendicularly oriented 
LDOS maps*i. Part of the modulation wavevectors ob- 
served in the experiments^ are also dispersive or energy 
dependent, this behavior should be related to the scat- 
tering of quasiparticles from defects including the ran- 
domly distributed stripes which break the translational 
invarianc o 23 i 24 . Here we would like to emphasize that 
our quasiparticles are excitations from the SDW/CDW 
(stripes) + dSC ground state. The information of the 
stripes is nonpertubatively included in wave functions of 
the quasiparticles. After the Hamiltonian is diagolized, 
there exists no extra term which couples the quasipar- 
ticle with the stripes. It should also be interesting to 
calculate the energy-dependent LDOS in the presence of 
stripes and to sec its difference as compared with the case 
for a pure dSC. In Fig. 2, we show the energy dependence 
of the LDOS at one of the hole (electron) accumulated 
stripes, namely at one of the minima (maxima) in Fig. 
1 (c) . It appears that there exist two small bumps within 
the superconducting coherence peaks, a signature of the 
presence of weak SDW and CDW orders. At the differ- 
ent sites, the LDOS also can be shown to have the similar 
features. 



FIG. 3: The spatial distribution of (a) the dSC order param- 
eter, (b) the staggered magnetization ,and (c) the electron 
density variation with a single strong impurity at site (0,0). 
(d) The LDOS maps at the energy E = 0.04. The param- 
eter values are the same as the Fig. 1, except the impurity 
scattering potential Vb = 100. 



Next let us introduce a nonmagnetic strong impurity 
like Zn into our system at site (0, 0) and investigate its 
effect. Here the on-site impurity potential is taken to be 
Vb = 100. From Fig. 3(a), one can easily see that the 
dSC order is suppressed and recovers to its bulk value at 
a length scale of £o, the superconducting coherence length 
(~ 5a), away from the impurity. Beyond this range, the 
weak stripe-modulation in the dSC order parameter still 
remains. In Fig. 3(b) we show that the impurity induced 
staggcrcd-momcnt of the SDW is zero at the impurity site 
and reach the maximum value at its four NN sites. The 
net induced local moment by the impurity corresponds 
to a local spin with S z =l/2 when U becomes stronger—. 
From Fig. 3(b), the SDW is clearly pinned at the impu- 
rity site with one of its ridges. The spatial profile of the 
electron density change is presented in Fig.3(c), and it 
exhibits a Freidcl-likc oscillation around the impurity. It 
is easy to see that the electron number reaches to zero or 
Srii = —0.85 at the impurity site, which is much larger 
than the amplitude ( 0.001) of the CDW stripes. This 
is the reason why the stripes in Fig. 1(c) are too faint 
to see in Fig. 3(c). Fig. 3(d) displays the LDOS map at 
energy E = 0.04. The intensity is zero at the impurity 
site and there are four peaks appear on the four NN sites 
around the impurity. The intensities of the two peaks 
at the site (0, ±1) are less than those at the site (±1, 0). 
The stripe structure in Fig. 1(d) still prevails, but it is 
again too weak to see here. 

It is well known that a sharp single resonance peak 
at zero bias^i in the LDOS at the NN sites of a uni- 
tary impurity appears for a pure dSC. In the presence of 
weak stripes, the LDOS as a function of energy at the 
NN site, (0,1) and (1,0) around the impurity are respec- 
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FIG. 4: The LDOS as a funtion of energy at the NN site (a) 
(0, 1) and (b) (1, 0) around the single strong impurity located 
at the site (0,0). The parameter values are the same as the 
Fig. 1, and the impurity scattering potential Vb = 100. 



tively displaced in Fig. 4(a) and (b). The zero bias peak 
previously obtained for U = is dramatically depressed 
by U 7^ and splits into two distinct peaks. This is be- 
cause stronger local AF order is induced near the strong 
impurity site and a larger SDW gap opens up locally that 
suppresses the LDOS close to the impurity. This makes 
the LDOS in the present case very different from that 
of a pure d-wave (or U — 0) case. The oscillation in 
the LDOS at negative bias below the left peak is orig- 
inated from the energy-dependent modulations induced 
by the impurity, and the large amplitude seems to come 



from the size effect. When U/g is not large enough to 
generate extended SDW, the weak and local AF order 
induced by the impurity may force the resonance peak 
just to split a litt lo if 26 . The shapes of split peaks at site 
(0, 1) and (1,0) are slightly different. This indicates that 
the fourfold symmetry for a pure dSC changes to twofold 
symmetry when the stripe phase is in presence. If the 
experimentally observerU^^ stripe structure is of mag- 
netic origin, wc predict that the LDOS should exhibit a 
double-peak feature near the zero bias. 

In conclusion, wc have studied a cuprate superconduc- 
tor with weak stripe- like modulations in the dSC, SDW 
and CDW order parameters without and with a strong 
nonmagnetic impurity. In the absence of the impurity, 
the LDOS exhibits two small bumps within the supercon- 
ducting coherence peaks, a signature of the presence of 
the competing AF order. The LDOS maps displays the 
same stripe modulation as the CDW along y-direction 
with periodicity 4a which is nondispersive or energy in- 
dependent. The components of the modulation wavevec- 
tors observed by the experiments which are energy de- 
pendent or dispersive should be attributed to the effect 
due to scatterings of quasiparticles from defect o 23 i 24 . In 
the presence of a strong impurity, wc predict that the 
LDOS at the NN sites of the impurity are strongly sup- 
pressed by the AF order and reveal a double peak struc- 
ture. Hopefully our theoretical results could be useful 
to future STM experiments performed on samples with 
weak stripe-like structures. 
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